Impurity state in the vortex core of d-wave superconductors: Anderson impurity 

model versus unitary impurity model 
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Using an extended Anderson/Kondo impurity model to describe the magnetic moments around 
an impurity doped in high-T c d-wave cuprates and in the framework of the slave-boson meanfield 
approach, we study numerically the impurity state in the vortex core by exact diagonalization of 
the well-established Bogoliubov-de Gennes equations. The low-energy impurity state is found to be 
in good agreement with scanning tunneling microscopy observation. After pinning a vortex on the 
impurity site, we compare the unitary impurity model with the extended Anderson impurity model 
■ by examining the effect of the magnetic field on the impurity state. We find that the impurity 

resonance in the unitary impurity model is strongly suppressed by the vortex; while it is insensitive 
^ \ to the field in the extended Anderson impurity model. 

. PACS numbers: 74.25. Jb, 72.15.Qm, 74.60.Ec, 73.20Hb 

Effects of nonmagnetic impurities, such as Zinc, doped in the Cu02 plane of high-T c materials have attracted 
much attention in recent years as it may help us to better understand the underlying mechanism of high-T c 
superconductijotyEI. The low-energy (near zero) impurity resonant state in d-wave superconductors was first pre- 
dicted in Refs-Ba based on a unitary scattering-potential impurity model and later successfully observed in a series of 
beautiful atomic-scale scanning tunneling microscopy (STM) experimentsETLl in the vicinity of individual Zn ions in 
Bi2Sr2Ca(Cui_ 2; Zn 2 ;)208+5. Since Zn ++ has no spin itself (therefore nonmagnetic), it is natural to treat the impurity 
as a point-like scalag potential scatterer of conduction electrons. Indeed, both energy position and four-fold symmetric 
spatial distributionEl of local density of states (LDOS) of quasiparticle can be explained consistently by theoretical 
calculations based on the i-matrix scattering potential theory in the unitary limito. Nevertheless, the sign and mag- 
i__T nitudc of the scattering potential as well as the particle- hole symmetry play a crucial role in determining the exact 
energy level of impurity resonant state relative to Eermi energy. The scattering potential chosen in various theoretical 
studies differ largely from —4.4 eV to 18.9 eVOllJ. The continuum ipinatrix theory, which assumed paxticle-hole 
symmetry and repulsive potential, predicted impurity state with energyB consistent with STM observation^ that the 
£^ ■ resonance lies slightly below the Fermi fevel; on the other hand, the particle-hole symmetry is not applicable to high-T c 
iy) | cuprates and theoretical calculationgj'Ej show that in the absence of suchpSymmetry only strong attractive potential 
. can produce impurity state with negative energy while repulsive potential!!!] gives positive energy contrary to exper- 
iment results. Furthermore, the unitary impurity model has difficulty in accounting for spectra distribution pattern 
of differential conductance at the resonance energy. The experiment shows that the zero bias peak is strongest on 
the Zn ion with local maximum on its second nearest neighbor Cu sitesu while the unitary impurity model concluded 
+^ \ that the spectral, weig ht on Zn atom is negligibly small and local maximum peaks are on nearest neighbor Cu sites 

S of the impurityolilniS. In order to reconcile the apparent discrepancy, the effect of the interlayer tunneling matrix 
elements, such as a blocking effect of BiO layerLH and a fork modelEj have been employed to reproduce the distribution 
i of spectral weight . At present, whether the spatial pattern of the tunneling conductance of the impurity resonance 
is intrinsic in the Cu02 plane or interlayer tunneling effect must be considered is still an open question for b 
oretical investigation and further experimental examination. Recently, several alternative theoretical modelsl 
have been suggested, which take the effect of the Kondo screening of the local magnetic moments induced around the 
Zinc atom into accountyjBostly motivated by the observation of nuclear magnetic resonance (NMR)E-3~li3 and neutron 
scattering experimentsEZHl3 that staggered magnetic moments are developed around Zn ion with total net spin 1/2. 
When Kondo spin dynamics of such moment is considered, by including an exchange interaction between the induced 
spin and the spin of conduction electrons in the model Hamiltonian, Polkovnikov et au produce an impurity state 
whose energy dependence and spatial pattern may fit well with STM spectra, without a strong scattering potential 
and the specific characteristic of interlayer tunneling matrix elements. 

Motivated by the above controversy between the unitary impurity model and the model invoking induced magnetic 
moments, in this work we intend to compare these two models in the case that a single vortex line is pinned by the 
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impurity. To describe the Kondo effect of the magnetic moments and the effect of vortex, we introduce an extended 
Anderson impurity model with a phase factor being dependent on site and vortex. In the absence of vortex, both 
models with appropriate parameters can generate low-lying impurity state. However, when a vortex is present, for the 
unitary impurity model, the pronounced LDOS peak generated around the Zn impurity is largely decreased or even 
destroyed, indicating the vanishing of the impurity resonance by the vortex core; while for the extended Anderson 
impurity model, the existence of the vortex center has a weak effect on the impurity state. Such a difference may be 
examined readily by STM experiments on B^Si^Ca^ui-zZn^Os+a samples subject to a strong magnetic field. 

In this work, we adopt an extended Hubbard model on a two dimensional (2D) lattice with nearest-neighbor (NN) 
hopping and NN pairing interaction to model the d-wave high-T c cuprates and the extended Anderson impurity model 
to describe the magnetic moments around the impurity (Zinc ion) doped in the CuC>2 plane. The model Hamiltonian 
is expressed as 



where 



H — Hdsc + H imp (1) 

H dsc = - XI *« c L c jV + ( Ai J c !t c ]i + h - c -) ~ A*clc i(T , (2) 

and 

H-wp = H mag + Vicl a c 0(T . (3) 

Here Hd sc is the BCS-like Hamiltonian of the host <i-wave superconductor described on a 2D lattice. refers 
to the NN sites and tij the hopping integral between site i and j. /i is the chemical potential. Aj,j is the bond 
paring potential defined as Ay = — V^CijCjj) with Vd the effective pairing strength between electrons on NN sites. 
iAmp represents the impurity Hamiltonian which includes both an on-site scattering potential term represented-by Vj 
and the term i? mag describing the local moments around the impurity. As suggested by NMR experimentsEJ~t3, the 
induced magnetic moments with net spin 1/2 are mainly located at the fcoir NN Cu sites of the central Zn impurity 
site Tq — (0, 0) while negligible right on the Zn site; therefore as in Refit], we assume that the magnetic impurity 
with an effective 1/2 spin may reside only in the NN sites and couples with conduction electrons on the NN sites. 
Namely, the H mag may be modeled by the extended Anderson impurity model with strong Hubbard repulsion Ud as 

ffmag = J2a e d d id a + U d d\d^d\di + J2 ri =T,a V h W ( e0i c \a d <j + h.c.) with t = ±x,±y unit vectors, where e d is the 
d-electron energy level and Vh is the hybridization of the moments with the conduction electrons. As a modification 
to the coupling term in the ordinary Anderson impurity model, we here introduce a site-dependent phase factor fa, 
which is found to be quite crucial in determining the energy level and even the existence of the impurity state if Vj is 
not so strong. Note that under the meanfield decoupling scheme of Ref.E3, the present model for iJ mag is essentially 
equivalent to the magnetic impurity(Kondo) model in Ref.B, where it is written as J^. KiS ■ s(i) with the summation 
being over the NN sites of the impurity, and S = Ylad da< J apdp/2 and s = c a a af3 c 0/^ representing respectively 
the spin operators of the magnetic impurity and the conduction electrons. This is because that this spin-exchange 
term can be transformed by the mean- field decoupling scheme in RefM KJ2i S ■ s(i) -> -3K/8J2Ax*4<7 d <? + h - c -}> 
where a complex Hubbard-Stratonovich field Xi = ( c t a da) is introduced as an effective hopping amplitude between the 
d-level and conduction electrons. At this stage, one can see clearly that the phase factor fa introduced in our model 
and that of the field Xi are closely related. Assuming the four- fold rotational symmetry with respect to the impurity 
site, there are four inequivalent arrangements of </>x, that is 



M, (4) 



with m = 0, 1,2,3, 0(±x) = (tt/2) =p (tt/2), and 0(±y) = ±(n/2). In Ref.EI, a nontrivial c£-wave pattern of Xi nas 
been found to have the lowest free energy saddle point, whose phase corresponds to m = 2, matching the underlying 
symmetry of local bond pairing potential Aq^. When the vortex exists, the symmetry of the underlying pairing 
potential would change from the local d-wave pattern to that responding to the winding of phase. Therefore, it 
appears reasonable to expect that <p might change from m = 2tom = lasa vortex is pinned at the impurity site. 
We find that such a response of <j> to the magnetic field is quite crucial in keeping the impurity state in the vortex 
core. 

In the treatment of the extended Anderson impurity model at low temperatures, Ud is assumed to be infinite, aa 
usual, which forbids double-occupancy of electrons on the d-level. Therefore, the slave-boson mean-field theorycj'EiJ 
can be applied as in Refs.uc3 where the d-electron operator is written as dJ a = f£b with f a the spin-carrying fermion 
operator and b the holon operator. Furthermore, the single occupancy constraint J2a f^f<r + °' f b = 1 should be obeyed. 
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At the mean-field level, the holon operators b and are approximated by a c-number 60 an d the constraint is enforced 
on average by introducing a Lagrange multiplier Aq. Accordingly, the mean-field -ff; mp becomes 



H imp = J2^dftfa+ ]T V h (i)(e^clf a + h.c.) + V lC l a coa + \v(bt-l) (5) 

with renormalized parameters = + Ao and Vh — Vft&o- By applying the self-consistent mean- field approximation 
and performing the Bogoliubov transformation, diagonalization of the Hamiltonian can be achieved by solving the 
following Bogoliubov-de Gennes(BdG) equations: 

where u n , v n are the Bogoliubov quasiparticle amplitudes with corresponding eigenvalue E n and Hij = —S(ri + r — 

T i)\ t i^id - ($i,id + < *j,»J v fc eV M ~ fityiddlM ~ S ( r i ~ r o)Vi] + ^i,i d ^d- Here, i,j represent the index of the 2D lattice 
sites with id as the index for the magnetic moment residing on the NN sites of the doped impurity. Ajj- is defined 
between a pair of NN sites on the 2D lattice and is calculated according to the self-consistent condition: 

n 

and 

bl = 1 - 2^{|<J 2 /(S„) + |<| 2 [1 - /(£„)]} (8) 

n 

Once the BdG Eq. ^is solved self-consistently, the quasiparticle spectrum can be obtained and the LDOS proportional 
to the differential tunneling conductance observed in STM experiments is given by 

piv^E) =-Y}K\ 2 f(E n -E) + \v?\ 2 f(E n +E)] (9) 

n 

where f(E) represents the Fermi distribution function. 

Unitary Impurity Model. — First, we study the variation of the impurity state when a vortex line is sitting on the 
impurity site within a unitary impurity model. Exact diagonalization is applied and the BdG equations are solved 
self-consistently. In studying the electronic structure of vortex lattice in d-wave-superconductors, a magnetic unit cell 
which accommodates two vortices is usually employed in numerical studiesE3~E3. Such a method has been applied to 
treat impurity effect in the mixed state of s-waveE3and d-wavet£l superconductors. The primitive translation vectors 
of the magnetic unit cell are R x = aN x x and R y = aN y y, where a is the lattice constant and will be set as unity. The 
pairing potential winds by An around the magnetic unit cell. The quasiparticle amplitudes ti" and v n are classified 
by the magnetic Bloch quasi-momentum k. 

The parameters we choose are fi = —0.2t, Vd — 2.2£ which give rise to = 0.274i (accordingly A max ~ t) and 
T c = 0.45i. We here intentionally take a relatively large amplitude of the energy gap A max with respect to real 
materials, because a large energy gap can sufficiently lower the LDDS peak corresponding to the vortex core states 
and in our case a broad structure is achieved near the vortex centered, which enables us to distinguish the LDOS peak 
of the unitary impurity from that corresponding to the vortex core states. The on-site attractive scattering potential is 
chosen as Vj = — lOt. In Fig. [j]the LDOS at the NN site of Zn site is plotted. In the absence of magnetic field(dashed 
line), a sharp peak is found at Eo/A. max — —0.02 corresponding to the impurity resonance. In the inset, the spatial 
distribution of the LDOS at resonance energy, p(r,Eo) is shown, indicating that the spectral weight is concentrated 
at the NN sites of the impurity while vanishingly small at the impurity site. These results are qualitatively consistent 
with previous numerical investigations. When a vortex is pinned right on the impurity site, from Fig. |l|(solid line) 
we can see that the LDOS peak at the NN site of Zn is significanly suppressed and mixes with the vortex core states 
without identification of impurity resonance state any more. As a strong scattering potential, the unitary impurity 
can drive the spectral weight on it mainly to its NN sites no matter vortex is present or not. However, in the absence 
of vortex, the LDOS at a specific energy is raised, leading to a resonance state; while in the presence of vortex, the 
LDOS on NN sites is increased in a wide energy range as seen in Fig. |l| and thus it is more reasonable to state that 
the vortex core states are strengthened by the impurity scattering. Broadening of the impurity resonance state by the 
magnetic field exists when the supercurrent of the vortex disturbs the energy spectrum of <i-wave superconductors via 
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Doppler shiftEI. However, in our case the smearing of the impurity resonance state is not induced by the supercurrent 
because the impurity is on the vortex center where supercurrent is zero. We attribute it to the nature of topological 
singularity of the vortex center where the phase of pairing potential varies arbitrarily, which makes the vortex center 
itself a strong scatter of electrons^. Therefore, it is the vortex that modifies the background electronic structure 
which is responsible for the impurity resonance and accordingly makes the impurity resonance indiscernible. 

Extended Anderson Impurity Model. — We now address the impurity state in the presence of Kondo screening 
of the magnetic moments described by the extended Anderson impurity model. First we consider the. case in the 
absence of applied magnetic field. We set m = 2 in Eq. (Q) in consideration of free energy minimizationO. The exact 
diagonalization is performed on a 24 x 24 square lattice, which is seen to be sufficiently large due to the localized 
feature of the impurity state. The scattering potential Vj is simplified as zero (The role of nonzero Vj will be discussed 
later). The bare d-level is set as = — At and the coupling strength Vh — 2i(Kondo screening regime). The value 
of the Lagrange multiplier Ao is determined by minimizing the free energy of the system as indicated in the inset 
of Fig. H and the corresponding value of b is calculated self-consistently according to Eq. (||). In the inset of Fig. 
||, AF as a function of Ao is plotted, where AF represents the change of the free energy when the local magnetic 
moments couple with conduction electrons. We find that for the chosen parameters the lowest free energy occurs 
when Ao — 5.6i with bo — 0.64. Figure || gives LDOS's on the Zn site, its NN site and NNN site. There is a sharp 
peak at the energy Eq/ A max = —0.04, which is well consistent with STM measurement. Right on the Zn site, the 
peak is strongest and local maximum peaks are on the NNN sites; while on the NN sites peaks are rather weak and 
at — Eq/ A ma „ — 0.04. The spatial distribution of LDOS at E Q is shown in the inset of Fig. ||, also in agreement with 
STM spectraE and a theoretical study based on another approach on an essentially equivalent modeB. Our analysis 
shows that the impurity state induced by the magnetic moment is a bound state with a very short attenuation length 

Recent STM measurements!^! of impurity states induced by Zn ++ and Mg ++ doped in YBa2Cu06.9 reported that 
the resonance energies are at (—10 ± 2) and (4 ± 2) meV, being not necessarily near zero bias. Based on our model, 
with the fixed coupling strength e.g. Vh = It, both the bare d-level e<j and on-site scattering potential Vi could adjust 
the resonance position of the impurity state. When Vj = and td ranges from to —At, Eo/A max varies from —0.18 
to —0.04. An attractive(negative) Vj drives the energy position below the Fermi level further (e.g. when V/ = —It, 
Eo/A max varies from —0.42 to —0.26 as ranges from to — At), while repulsive one makes the resonance energy 
approach to or even across the Fermi level and become positive (e.g. when Vi — t, Eq/ A rnax — 0.25 with = =-4i). 
However, too large Vi (S> t) can affect the spin-induced impurity state strongly and double resonance peaksEj are 
found in the unitary limit, which disagrees with the STM spectra. 

When a vortex is pinned on the impurity site, we expect that there is a rearrangement of </>^ according to the 
winding of the phase of the order parameter so that the phases of the Hubbard-Stratonovich field Xi an( i the local 
bond pairing potential between NN sites Ao,i have the best match. Therefore, m will vary from 2 to 1 when a vortex is 
present. This is quite important as other choices of m (0, 2, 3) are found to result in the disappearing of the impurity 
state. The numerical results are shown in Fig. [|, which plots the LDOS as a function of energy on the impurity 
site. A strong single peak is seen near the Fermi level, indicating that the impurity state still survives under the field. 
Being different from the case without vortex, the peak is lowered slightly and moved to positive bias Eq/ A max = 0.02, 
while the spatial pattern of the impurity state is not modified(see the inset of Fig. || ). 

In summary, we have compared the extended Anderson impurity model with the unitary impurity model by studying 
numerically the impurity state of a d-wave superconductor as a vortex is pinned on the impurity site. We find that the 
impurity resonance governed by the scattering potential mechanism is sensitive to the presence of vortex, while the 
impurity bound state generated from the magnetic moments mechanism seems rather robust to the field. Therefore, 
it may be helpful to clarify the dominant mechanism of the impurity state in Zn or Mg impurities doped cupates by 
examining the existence of the impurity resonance, the variation of its energy, and even the STM spatial pattern in 
the presence of a strong magnetic field. 
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FIG. 1. LDOS versus energy at an NN site within the unitary impurity model for the cases (i) without vortex(dashed 
line), (ii) a vortex pinned at the impurity site (solid line), and (3) a vortex core without impurity (dotted line). The bulk 
density of states is also plotted for reference (dash-dotted line). Inset shows the spatial distribution of the LDOS at energy 
Eo/Amax = —0.02 in a 16 x 16 region for case (i). 



FIG. 2. LDOS versus energy at the impurity site (solid line), NNN site (dashed line), and NN site (dotted line) under 
zero field within the extended Anderson impurity model. The bulk density of states is also plotted for reference (dash-dotted 
line). Inset (a) shows the variation of AF as a function of Ao and (b) gives the spatial distribution of the LDOS at energy 
E /A max = -0.04. 



FIG. 3. Same as Fig. |2] except that a vortex is pinned on the impurity site. Inset (b) gives the spatial distribution of the 
LDOS at energy E /A max = 0.02. 
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